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A  mathematical  model  of  the  Sea  King  Mk  50  helicopter,  as  used  in  the  Anti-Submarine 
Warfare  (ASW)  role,  has  been  developed  at  ARL  to  run  on  the  Elxsi  6400  computer.  To 
validate  this  model,  extensive  flight  trials  have  been  conducted  by  the  RAN.  This 
document  provides  a  catalogue  of  the  many  flight  trials  data  files,  shows  how  to  access  and 
process  the  flight  data,  and  then  how  to  run  the  mathematical  model  with  inputs  obtained 
from  the  flight  data. 
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1 .  INTRODUCTION 

A  mathematical  model  of  the  Sea  King  Mk  50  helicopter,  as  used  in  the  ASW  (Anti- 
Submarine  Warfare)  role,  has  been  developed  by  ARL  to  a  Royal  Australian  Navy  (RAN) 
task  requirement.  This  model,  which  was  originally  developed  on  a  DEC  System  10 
computer  using  the  simulation  language  "CSMP-10(ARL)"  (Refs  1-3),  has  been  described 
in  general  terms  in  Refs  4  and  5.  Full  descriptions  of  the  main  components,  namely  the 
Aerodynamics/  Kinematics,  Control  Systems,  and  Cable/Sonar  may  be  found  in  Refs  6-11. 
A  large  data  bank  derived  from  an  extensive  series  of  Sea  King  flight  trials  (Ref.  12)  is 
presently  available,  the  processing  of  which,  together  with  the  application  of  filtering  and 
kinematic  consistency  checking  procedures,  has  been  described  in  Refs  13-16.  The  main 
purpose  of  the  trials  was  the  acquisition  of  data  for  the  validation  of  the  mathematical  model 
(Ref.  17). 

Assuming  a  basic  knowledge  in  using  the  Elxsi  6400  computer,  the  use  of  the  computer 
programs  associated  with  the  Sea  King  model  is  described  in  three  user  manuals.  Part  I 
(Ref.  18  )  shows  how  to  set  up  the  model  and  run  it  in  its  basic  modes  without  dunking 
sonar.  The  model  is  first  run  in  ASW  mode  as  a  means  of  trimming  the  aircraft,  i.e.  'flying 
to  trim',  and  then  in  either  ASW,  ASE  (Auto  Stabilizing  Equipment),  or  pilot  modes  to 
simulate  a  desired  manoeuvre.  Part  II  (this  one)  provides  a  catalogue  of  the  many  flight 
trials  data  files  (not  included  in  Ref.  12),  shows  how  to  access  and  process  the  flight  data, 
and  then  how  to  run  the  mathematical  model  with  inputs  obtained  from  the  flight  data.  Part 
III  (Ref.  19)  shows  how  to  use  the  dunking  sonar  model  and  demonstrates  the  use  of  a 
cable  graphics  program. 

2 .  RETRIEVAL  OF  FLIGHT  DATA  FROM  MAGNETIC  TAPE 

The  unprocessed  flight  trials  data  files  are  available  on  ARL  magnetic  tape  Ml 35.  An 
example  is  given  below  of  how  to  obtain  these  files,  a  catalogue  of  which  is  given  in 
Appendix  A. 

To  retrieve  these  files,  a  request  is  first  made  to  the  operator  for  tape  Ml 35  to  be 
mounted:1 


iMOOHTTAPK  M135  -W 


***Fran  operator  at  11:32:  Mag  tape  mounted  on  tapel 

:  TAPES 

Device  Status 

devioe  user  volume  type 

tapel 

tape2  ML35  ANSI 

tape3 

No  outstanding  mount  requests  for  user  ae.amey 


density  ring 


6250fcpi  No 


1  For  computer  terminal  input  included  in  (his  document,  messages  typed  by  the  user  are  shown  in  bold 
type. 


1 


In  order  to  copy  one  or  more  of  these  files  from  tape  Ml 35  to  the  current  directory, 
shellfile  REDATA  is  available  on  tape  M228  (see  Appendix  A  of  Ref.  18).  A  listing  of  this 
file  is  shown  below: 

:LIST  REDATA 

- SHEUFIIE  TO  RESTORE  LIST  OF  FIIES  FRCM  MTAPE 

parm  file  +list  +req 

- Init  shellvariables, 'files '“pathname, 'count'  is  a  counter 

set  filelist  ' '  ^declare 
set  count  1  -Kteclare 

set  filelist  teat  /subdotains/arl/user/ae.williams/sealdjng/flight_data/  & 

[file  [count]]]  +append 

-  Start  of  loop  to  add  list  of  files  to  pathname 

label  loop 

set  count  [eval  count+1] 

-  -  If  another  file  is  there  add  to  pathname  otherwise  restore  files 
if  [file  [count]]  then 

set  filelist  [cat  ",/subdCfnains/arl/user/ae.williams/seaking/flic£it_data/"  & 
[file  [count]]]  +append 

goto  loop 
else 

goto  restore 
end  if 

—  Restore  files 
label  restore 

restore  [filelist; noquote]  vol“M135  merge=f lat  SEQ=1  -uni  4cre 

Below  is  an  example  showing  how  REDATA  is  used  to  restore  flight  data  files  15018 
and  15019,  assuming  the  user  has  mounted  tape  Ml 35. 

.-REDATA  15018,15019 

15018 

15019 

**********  RESTORE  SIM4ARY  ********** 

2  file3  restored. 

0  files  not  restored. 

0  directories  restored. 

0  directories  not  restored. 

All  requests  were  found  on  tape. 


3 .  PROCESSING  THE  FLIGHT  DATA 

The  raw  flight  data  is  processed  using  the  data  processing  program  REFINE,  (see  Refs 
13-16).  The  program  requires  the  calibration  file  CAL86  (see  Appendix  B),  which,  as  well 
as  containing  details  of  the  instrument  calibrations,  also  contains  labels,  plot  limits,  and 
characteristics  of  various  filters  used.  Since  REFINE  was  originally  documented  in  Ref. 
13,  a  digital  notch  filter  has  been  added  (Ref.  16).  However,  the  processing  time  is 
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increased  dramatically,  and  it  is  usually  only  used  for  special  cases.  Both  REFINE  and 
CAL86  are  available  on  ARL  magnetic  tape  M228. 

For  the  flight  data  file  15018,  selected  here  as  an  example,  REFINE  is  run  as  follows: 


:CCFZIC 

INPOT  DATA  FILQMC  -  15018 
".DAT"  COTPOT  FILENAME  (w/o  ext)  « 

TITLE  (2  lines  of  60  chrs) 

:  1501 8  -  Bow  -  Pedal  Pulse  poet  -  ASS  on  -  Fit  4 


:  Is  daisy,  5s  liadt 

ARE  ASSIGNED  BLK  NUMBERS  REQRD  :  T 

CALIBRATION  FTLEJWC  -  CALB6 

FLIGHT  NLMER  -  4 

DCES  CHANNEL  18  t€ASUR£  TCBQUE  :  N 

IS  PRE-PROCESSING  WITH  A  NOTCH  FILTER  RE QRD  :  N 

OUTPUT  INTERVAL  (in  60'ths  of  sec;  e.g.  12  for  0.2  sec)  •=  1 

STARTING  TOE  DEIAY,  NO.  BUG  FIRST  IGNORED  -  1 

TOC  LIMIT  »  5 

IS  FILTERING  REQRD  :  T 

IS  OOTPOT  CF  FILTER  CHARACTERISTICS  REQRD  :  V 
ARE  DIGITAL  FILTER  DELAY  ADJUSTMENTS  RECBD  :  Y 
IS  3CC0HING  RECPD  :  N 

ARE  INSTRIMNr  «  ANALOGUE  FILTER  DELAY  ADJUSTMENTS  REQRD  ;  Y 

ARE  SCAIFS  AND  OFFSETS  REQRD  :  Y 

APE  PLOT  LIMITS  FEQFD  :  T 

ARE  DROP-COTS  TO  BE  COFRECTED  ;  Y 

APE  ALL  CHAMELS  RECFD  :  II 

OM8CL  NUMBERS  RECK)  FOR  FGUCWING  GROUPS 

(Set  first  value  -we  If  excluding) 

[Set  first  value  to  -100  if  none  excluded) 


Varies  from  1  to  6 

Ch.  18  measures  torque  or  yaw  attitude 
Used  in  special  cases 

Found  with  experience 
Found  with  experience 


Instrumentation  data 
Bocm  Calculations 
Blade  angles 
Euler  angles 


<  1  to  33) 
(34  to  40) 
(41  to  45) 
(46  to  47) 


Kinematic  consistency  (48  to  70) 
DATE  -  O9-Jan-07 
TIME  -  14:56:57 


1.2.  3.4,8,  9,11,12,13.14 ,15.16,17 ,21.  22.  27.  29,  30 

41,42,43,44,45 

46,47 

54,55,56,57,58,59 


IS  KINEMATIC  CONSISTENCY  ITERATION  INFORMATION  REQRD  :  N 

ARE  ANY  FAULTY  DATA  TO  BE  REPIACH)  BY  ALTERNATIVELY  DERIVED  DATA: 

No.  tine  correct  ions  -  0 

No.  blocks  replaced  *  0 

No.  drop-out  corrections  »  265 


Fortran  program  executed  STOP  statement  1 
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The  above  run  creates  the  files  15018.ERR  and  15018.DAT.  File  15018.ERR  contains 
a  list  of  errors  encountered,  including  time  corrections,  blocks  replaced,  and  drop-outs. 
File  15018.DAT  contains  the  processed  data  in  a  format  suitable  for  use  as  input  to  the  post 
processing  program  TRANS  (see  Sections  4.2  and  5). 

It  should  be  noted  that  suspect  angular  measurements  were  recorded  in  some  flights, 
particularly  of  yaw  angle  and  yaw  rate.  These  can  be  detected  by  comparing  with  values 
derived  through  kinematic  consistency  checking  (Ref.  15). 1  The  faulty  measurements  may 
then  be  replaced  by  the  alternatively  derived  (by  integration  or  differentiation  of  other 
measured  quantities)  values  (see  Appendix  B  of  Ref.  15). 

4 .  RUNNING  THE  SEA  KING  MODEL  WITH  FLIGHT  DATA  INPUTS 

4.1  Standard  Input  Files 

As  described  in  Section  5.1  of  Ref.  18,  whenever  the  model  is  run,  the  following  three 
files  are  required  as  input: 

BOMMP.IN  -  Non-interactive  command  file  for  BOMMP 

DATA.HEL  -  Helicopter  input  data,  mainly  in  NAMELIST  form 

?????.MOD  -  Helicopter  model  information  in  form  of  configuration, 

parameter,  and  function  statements  -  must  have  5  character 
name  with  .MOD  extension 

Standard  helicopter  input  files  and  corresponding  model  files  for  ASW  mode  are 
available  at  All-Up-Weight  (AUW)  values  of  16600,  17800,  and  19200  lb.  Where 
conditions  need  to  be  represented  more  specifically,  the  input  data  file  (DATA.HEL)  can  be 
edited  and  the  model  file  can  be  created  by  re-trimming  in  ASW  mode  (see  Section  5.2  of 
Ref.  18). 

4.2  Flight  Control  Inputs 

4.2.1  Pedal  Inputs 

For  the  file  15018  processed  in  Section  3,  the  AUW  is  given  as  16640  lb  in  Appendix 
A.  Because  this  is  close  to  the  value  16600  lb,  the  standard  files  16600.HEL  (renamed  to 
DATA.HEL)  and  16ASW.MOD  are  used.  The  program  SKMODE  (see  Section  5.3  of 
Ref.  18)  is  first  run,  as  shown  below,  to  obtain  the  ASE  mode  file  16ASE.MOD: 


: SKMODE 

INPUT  FIIE  (ASW  MX®  .MOD  FILE)  :  16 ASW 
MODEL  FILE  REQUIRED  (ASW,ASE,PIL  OR  ALL)  :  ASE 
CUIPUT  FILE  (ASE  MX®  .MX)  FII£)  :  16ASE 

??ASE.MX>  TITIZ  :  SEA  KING  -  Hover  -  3000  ft  -  ASE  MODE  -  166001b  AUW 


1  Flight  data  files  which  have  previously  been  found  to  have  faulty  data  are  indicated  in  Appendix  A. 
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The  file  BOMMP.IN  is  then  edited,  replacing  19ASW  by  16ASE,  to  give: 


•■LIST  BOMMP.IN 

L0G2:16ASE_I 

CCN 

PAR 

FUN 

MAN 


Finally,  a  special  file  containing  the  flight  data  control  inputs  is  created  using  the 
program  TRANS  (Refs  2  and  3),  which  is  available  on  ARL  magnetic  tape  M228.  This  file 
is  read  in  by  the  Sea  King  model,  any  deviation  from  the  initial  control  position  being  added 
to  the  model  trim  condition.  Using  the  data  file  15018.DAT  (created  in  Section  3)  as  input 
to  TRANS,  the  file  15018.COL  (see  Appendix  C)  is  created  for  use  as  input  to  the  Sea 
King  model  as  shown  below: 

:  TRANS 

1TRPNS  version  date  ll-MAR-86] 

I/P  FII£NAME  -  15018 

15018  -  Hover  -  Pedal  Pulse  Port  -  ASE  cn  -  Fit  4 

I/P  FI1E  RECORDED  CN  09-dan-87  AT  14:56:57 
INTEEN  INT  a  -0000E+00;  RUN  CPU  TIME  =  24.98  SEC. 

TIME  FRCM  .0000E+00  TO  5.0000E+00  IN  STEPS  OF  I.6667E-02 

*PRC 

PRINTING  IN  OOUJFNS  : 

BIKS 

276, 277,72, 278,527,292, 293,295, 294 

IS  O/P  TO  Try  RECFD  :  H 
*GOE 

**  RUNNDJ3  ** 


*EXI 

Having  prepared  all  the  relevant  files  (i.e.  BOMMP.IN,  16ASE.MOD,  DATA.HEL, 
and  15018.COL),  the  model  may  now  be  run.  Depending  on  the  type  of  input,  a  number 
of  switches  need  to  be  turned  on  or  off.  The  example  below  shows  how  this  is  done  using 
the  PARameter  command,  for  the  case  here  (see  Appendix  A),  where  there  is  a  pedal  input 
with  the  Automatic  Stabilizing  Equipment  (ASE)  on: 
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:SEAKING86 

MAX  BLK  NO.  =  500 


MAX  NO.  OF:  I  S  T1  BIKS,  U  BLKS,  F  BIKS  =  100,25,25 

SEA  KING  -  HOVER  -  3000  ft  -  ASE  MODE  -  16600  lb  AUW 

•PAR 

PARAMETERS  : 

BLK,  PI,  P2,  P3 

57.1 

74.1 

79.1 

87.1 
90,0 
241,-1 

99.1 
123,0 
274,500 

•TNT 

INTH3J  P ARAMS;  IOWER,  UPPER,  INTERVAL  =  0,5,0.02 

•OUT 

O/P  BIKS 

276,277,72, 278, 292, 293, 295,  296 
282,283, 284,279, 280,281,287,288 
16,17,18,212,290 

O/P  PARAMS;  %  CHANGE  REQRD,  INTERVAL  =*  0.0001,0.1 

•LOG3: 15018_M 

MODEL  O/P  TO  1033:15018  -MOD 

•STO 

15018  .MOD  NOT  ON  DSK 


Read  cyclic  stick  position 
Read  yaw  controls 
Read  yaw  push  ■ pull-rod  position 
Read  collective  slick  position 
Do  not  read  blade  angles 
Do  not  read  blade  angles 
Display  monitor  output 
Turn  off  yaw  ASE  channel 
Set  yaw  ASE  channel  on  in  500  s 


CON,  PAR,  TON,  CR  ALL  :  A 
•LOG2 : 15018_I 

MODEL  I/P  FROM  LCG2: 15018.  MOD 

•LOG1 : 151 8A_0 

BIOX  O/P  TO  1031: 1518A.DAT 

•GOE 

1518A.DAT  NOT  ON  DSK 
**  RUNNING  ** 

1518A. HEL  NOT  ON  DSK 
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(s) 

(Ft) 

(Ft/m) 

(Knots) 

Time 

Alt 

RoC/D 

Speed  Slip 

%Torq 

.0 

3000. 

-.1 

1 

O 

1 

o 

o 

76 

%Colctv  SCyclic 

F-A  Lat 

-103.14  -19.82  -1.10 

A.S.E.  Channels  : 

PITCH  -  On 
ROLL  -  On 
YAW  -  Off 
ALT  HOLD  -  Off 

A.S.W.  Mode  : 

TRAN 

(Cyclic  Trim  -  RELS] 


.0 

3000. 

-.1 

-.0 

-.00 

76 

1.0 

3000. 

.0 

-.0 

.02 

76 

2.0 

3000. 

39.6 

-.2 

.62 

86 

3.0 

3001. 

81.7 

-.1 

1.23 

71 

4.0 

3003. 

111.6 

-.1 

1.90 

83 

5.0 

3005. 

150.1 

-.1 

2.73 

77 

RUN  CPU  TIME  :  15.36  Sec. 


»EXI 

Fortran  program  executed  STOP  statement  0 

The  above  example  creates  the  output  file  1518A.DAT,  which  can  then  be  used  as  input 
to  program  TRANS,  to  produce  suitable  output  in  tabular  or  graphical  form  (see  Section  5). 

Further  explanation  of  the  switches  set  using  the  PARameter  command  is  now  made. 
When  reading  flight  control  inputs,  the  cyclic  stick  and  collective  stick  positions  are  always 
read  (Blks  57  and  87),  as  well  as  a  yaw  channel  input  (Blk  74).  For  pedal  inputs,  the  yaw 
push-pull-rod  movements  are  read  (Blk  79  =  1),  whereas  for  other  inputs,  the  pedal 
position  is  read  (Blk  79  =  0).  This  is  because  the  yaw  control  system  is  not  modelled 
correctly,  and  variables  required  to  model  the  system  were  not  measured  (Ref.  17).  In  the 
case  of  a  pedal  input,  a  pedal  position  is  derived  from  the  push-pull-rod  position,  assuming 
no  input  from  the  ASE,  and  making  no  allowance  for  the  faulty  modelling.  For  'ASE  on' 
modelling,  the  yaw  channel  should  be  turned  off  (Blks  123  and  274),  since  ASE  inputs  are 
already  included  in  the  yaw  push-pull-rod  flight  data.  The  Sea  King  model  also  has  the 
facility  to  read  in  flight  data  blade  angles  (Blks  90  and  241),  instead  of  control  positions, 
thus  by-passing  the  modelling  of  the  flight  control  system. 

4.2.2  Non-Pedal  Inputs 

Below  is  an  example  of  how  the  model  is  normally  run  for  non-pedal  inputs,  file  15019 
representing  a  collective  input  manoeuvre: 
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:SEAKING86 

MAX  BIX  NO.  -  500 

MAX  MO.  CF:  I  S  T1  BLKS,  U  BLKS,  F  BLKS  =  100,25,25 

SEA  KING  -  HOVER  -  3000  ft  -  ASE  MODE  -  16600  lb  AUW 

•PAR 

PARAMETERS  : 

BIX,  PI,  P2,  F3 

57.1 

74.1 
79,0 

87.1 
90,0 
241,-1 

99.1 

123.1 
274,0 

*INT 

INTBGN  PARAMS;  BOWER,  UPPER,  INTERVAL  -  0,5,0.02 

•OUT 

O/P  BIXS 

276,277,72,278,292,293,295,296 

282,283,284,279,280,281,287,288 

16,17,18,212,290 

O/P  PARAMS;  %  CHANGE  RBQRD,  INTERVAL  =  0.0001,0.1 

•LOG3 : 15019_M 

MODEL  O/P  TO  IOG3.-15019.M3D 

•STO 

15019  .MOD  NOT  CN  DSK 


Read  cyclic  slick  position 

Read  yaw  controls 

Read  pedal  position 

Read  collective  slick  position 

Do  not  read  blade  angles 

Do  not  read  blade  angles 

Display  monitor  output 

Turn  on  yaw  ASE  channel 

Set  yaw  ASE  channel  on  at  time  =  Os 


cat,  PAR,  FUN,  OR  ALL  :  A 
•LOG2 : 15019_I 

MODEL  I/P  FROM  LOG2: 15019. MOD 

•LOG1 : 1519A_0 

BLOCK  O/P  TO  D0G1 : 1519A.DAT 


•GOE 

1519A.DAT  NOT  CN  DSK 
**  FUNNING  ** 
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1519A. HEL  NOT  CN  DSK 


(s)  (Ft)  (Ft/m)  (Knots) 

Time  Alt  RoC/D  Speed  Slip  %Torq 
.0  3000.  -.1  -.0  -.00  76. 

%Colctv  %Cyclic 

F-A  Lat 
-103.14  -19.82  -1.10 

A.S.E.  Channels  : 

PITCH  -  On 
ROLL  -  On 
YAW  -  cn 
ALT  HOLD  -  Off 

A.S.W.  Mxte  : 

TRAN 


[Cyclic  Trim  - 

.0  3000. 

RELS] 

-.1 

-.0 

-.00 

76 

1.0 

3000. 

146.1 

.1 

.02 

87  : 

2.0 

3000. 

394.2 

.4 

.25 

93 

3.0 

3001. 

544.8 

.6 

.96 

89 

4.0 

3003. 

653.1 

.8 

2.22 

90 

5.0 

3005. 

736.9 

1.0 

3.58 

90 

RUN  CPU  TIME  : 

15.14 

Sec. 

*EXI 

Fortran  program  executed  STOP  statement  0 


The  output  data  stored  as  file  1519A.DAT  may  then  be  processed  using  TRANS  (see 
Section  5). 

4.3  Engine  Cut 

During  the  Sea  King  flight  trials  (Ref.  12),  a  number  of  tests  were  done  simulating  an 
engine  failure,  where  one  engine  was  closed  down  to  idle,  while  the  other  engine 
automatically  increased  power.  In  order  to  simulate  these  tests  with  the  Sea  King  model,  a 
number  of  changes  had  to  be  made  to  a  model  file,  while  leaving  the  Fortran  source  files 
unchanged. 

For  the  file  19092  considered,  Appendix  A  gives  the  AUW  as  19100  lb,  the  aircraft 
forward  speed  as  80  KIAS,  and  a  cut  of  Engine  2.  Being  close  to  the  AUW  of  19200  lb, 
the  standard  files  19200.HEL  (renamed  to  DATA.HEL)  and  19ASW.MOD  are  used. 

On  inspecting  the  processed  flight  data  (in  which  Channel  18  measures  torque),  it  was 
found  that  the  aircraft  was  actually  at  a  forward  speed  of  82  kn  in  a  4  deg  climb.  The 
following  steps  are  required  to  achieve  a  trimmed  'ASE  mode'  model  file . 

Initially,  the  'ASW  mode'  is  used  to  'fly  to  trim'  with  a  forward  velocity  of  82  kn 
(138.4  ft/s)  in  level  flight.  This  is  demonstrated  in  the  example  below,  where  set  forward 
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speed  (Blk  51)  is  ramped  from  0  to  138.4  over  a  200  second  interval,  and  then  set  constant 
at  138.4  for  a  further  200  seconds  to  stabilize  the  trim: 

:Seaking86  Having  lirst  specified  19ASW.MOD  as  model  input  file  in  BOMMP.IN 

MAX  B1K  tC.  -  500 

MAX  NO.  CF:  I  S  T1  BLKS,  U  BLKS,  F  BLKS  =  100,25,25 

SEA  KING  -  HOVER  -  3000  ft  -  ASW  MODE  -  19200  lb  AUW 

*TIT 

TITLE  (LIMIT  60  CHRS) 

SEA  KING  -  82  ten.  Level  Flight  -  3000  ft  -  ASW  MODE  -  19200  lb  AUW 
*FUN 

FUNCTIONS  -. 

blk  no.  =  51  Set  aircraft  foiward  velocity  (ft/s) 

COORD  PAIRS  : 


0,0 

COORD  PAIR  (  -OOOOE+OO,  .OOOOE+OO)  DELETED 

200.138.4 

400.138.4 

BLK  NO.  = 

MODEL  OCRPLETE 

*INT 

INTK3J  PARAMS;  LOWER,  UPPER,  INTERVAL  =  0,400,0.02 

‘OUT 

O/P  BLKS 

2  Arbitrary  output  required  lor  intermediate  step 

O/P  PARAMS;  %  CHANGE  REQRD,  INTERVAL  “  0.0001,100 

*LOG3 :  TRIM2_M 

MODEL  O/P  TO  LOG3:TRIM2  .MOD 

*STO 

TRIM2.MCD  NOT  CN  DSK 


O0N,  PAR,  FUN,  OR  ALL  :  A 
*LOG2:TRIM2_I 

MODEL  I/P  FROM  LOG2 :  TRIM2 .  MOD 

*LOGl : TRIM2_0 

BLOCK  O/P  TO  LOG1: TRIM2.DAT 

*GOE 
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TRIM2  .DAT  NOT  CN  DSK 


»*  FUNNING  ** 


TRIM2.HEL 

NOT  CN  DSK 

(S)  (Ft) 

(Ft/m) 

(Knots) 

Tine  Alt 

FoC/D 

Speed 

Slip 

%Torq 

.0  3000. 

.3 

.0 

.00 

90 

%Colctv  %Cyclic 

F-A  Lat 


-94.31  -25.06  -2.49 

A.S.E.  Channels  : 

PITCH  -  On 
POLL  -  On 
YAW  -  On 
ALT  HOLD  -  RAD 

A.S.W.  Mode  : 

TRAN 


[Cyclic  Trim  ~~ 

.0  3000. 

ENGE] 

.3 

.0 

.00 

90. 

1.0 

3000. 

.3 

.0 

.00 

90. 

2.0 

3000. 

.3 

.1 

.00 

90. 

3.0 

3000. 

.3 

.2 

.02 

90. 

4.0 

3000. 

.1 

.4 

.04 

90. 

5.0 

3000. 

-.0 

.6 

.07 

90. 

398.0 

3000. 

.2 

81.9 

-.00 

45. 

399.0 

3000. 

.2 

81.9 

-.00 

45. 

400.0 

3000. 

.2 

81.9 

-.00 

45. 

RUN  CPU  TIME  : 

7  Min.  38. 

03  Sec. 

*RET 

*STO 

TRIM2.M0D  NOT  CN  DSK 
CON,  PAR,  FUN,  OR  ALL  :  A 
*EXI 

Fortran  program  executed  STOP  statement  0 

To  now  trim  the  model  in  a  4  deg  climb,  the  height  gain  over  a  200  second  interval  is 
first  calculated  (=  138.4  x  200  tan  4°  =  1935.6  ft/s).  The  aircraft  set  height  (Blk  80)  is  then 
ramped  from  its  starting  value  of  3000  ft  to  the  final  value  of  4935.6  ft  (=  3000  +  1935.6) 
over  this  200  second  interval,  while  keeping  aircraft  set  speed  (Blk  51)  constant  at  138.4 
ft/s: 
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:SEAKING86 

MAX  BIX  MO.  =  500 


Having  first  specified  TRIM2.MOD  as  model  input  file  in  BOMMP.tN 


MAX  NO.  CF:  I  S  T1  B1KS,  U  BLKS,  F  BLKS  =  100,25,25 

SEA  KING  -  82  ten.  Level  Flight  -  3000  ft  -  ASW  MODE  -  19200  lb  AUW 

*TIT 

TITIE  (LIMIT  60  CHRS) 

SEA  KING  -  82  ten,  4 deg  Climb  -  3000  ft  -  ASH  MODE  -  19200  lb  AUW 
‘FUN 

FUNCTIONS  : 

bix  no.  =  51  Set  aircraft  forward  velocity  (ft/s) 

COORD  PAIRS  : 


0,138.4 

COORD  PAIR  (  .0000E4OQ,  .OOOOE+OO)  DELETED 

bix  no.  “  80  Set  aircraft  height(ft) 

coord  pairs  : 

200,4935.6 

BIX  NO.  = 

MODEL  COMPLETE 

•TNT 

INTEEN  PARAMS;  LOWER,  UPPER,  INTERVAL  =  0,200,0.02 

‘OUT 

O/P  BLKS 

2  Arbitrary  output  required  tor  intermediate  step 

O/P  PARAEG;  %  CHANGE  REQRD,  INTERVAL  =  0.0001,100 

*LOG3:4DCLM_M 

MODEL  O/P  TO  LOG3:4DCLM.M3D 

*STO 

4DCIM.MCD  NOT  CN  DSK 


OON,  PAR,  FUN,  OR  ALL  :  A 
*LOG2 :  4DCLM_I 

MXIEL  I/P  FRCM  L0G2:4DCIM.MCD 

•LOG1 ;  4DCLM_0 

BLOCK  O/P  TO  L0G1: 4DOM.DAT 
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*GOE 

4DdM.DAT  NOT  ON  DSK 
**  RUNNING  ** 

4DCIM.HEL  NOT  ON  DSK 

(s)  (Ft)  (Ft/m)  (Knots) 

Tine  Alt  Ro C/D  Speed  Slip  %Torq 
.0  3000.  .2  81.9  -.00  45. 

%Colctv  ICyclic 

F-A  Lat 
-116.65  -4.10  -3.95 

A.S.E.  Channels  : 

PITCH  -  On 
POLL  -  On 
YAW  -  On 
ALT  HOLD  -  RAD 


A.S.w.  Mode  : 

TRAN 

(Cyclic  Trim  -  EDGE] 


.0 

3000. 

.2 

81.9 

\ 

© 

o 

45 

1.0 

3001. 

192.8 

82.0 

.51 

53 

2.0 

3006. 

394.4 

82.2 

-.73 

61 

3.0 

3014. 

492.8 

82.1 

-.90 

57 

4.0 

3022. 

515.5 

82.0 

-.40 

57 

5.0 

• 

• 

3031. 

499.8 

81.8 

-.02 

54 

• 

• 

198.0 

4872. 

580.8 

81.9 

.00 

60 

199.0 

4881. 

580.9 

81.9 

.00 

60 

200.0 

4891. 

580.8 

81.9 

-.00 

60 

RCJN  CPU  TIME  :  3  Min.  53.37  Sec. 

*RET 

*STO 

4DOM.MCD  NOT  CN  DSK 
OCN,  PAR,  FUN,  CR  ALL  :  A 

*EXI 

Fortran  program  executed  STOP  statement  0 

The  model  file  4DCLM.MOD  obtained  above  is  then  used  as  input  to  program 
SKMODE  in  order  to  obtain  the  ASE  model  file  4DASE  as  follows: 
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:SKMODE 

INPUT  FILE  (ASW  MCDE  .MX)  FILE)  :  4DCLM 
MCDEL  FILE  REQUIRED  (ASW,ASE,PIL  OR  ALL)  :  ASE 
OUTPUT  FILE  (ASE  MXE  .MX)  FILE)  :  4DASE 

??ASE.MX)  TITLE  :  SEA  KING  -  82kn,  4 deg  Climb  -  ASE  MODE  -  192001b 

This  model  file  can  then  be  used  as  input,  as  shown  below,  to  obtain  the  results  shown 
in  Appendix  D: 

:  seaking86  Having  first  specified  4 D ASE. MOD  as  model  input  file  in  BOMMP.IN 

MAX  BIX  to.  =  500 

MAX  NO.  CF:  I  S  T1  BIXS,  U  BLKS,  F  BLKS  =  100,25,25 

SEA  KING  -  82kn,  4deg  Climb  -  ASE  MX®  -  192001b 

*TIT 

TITLE  (LIMIT  60  CHRS) 

SEA  KING  -  82  kn,  4deg  Climb  -  Engine  Cut,  ASE  MODE  -  19200  lb 
*CON 

OCNFIGURATICNS  : 

BIX,  TYPE,  Bl,  B2,  B3 

316, +,-183, 326;  trq  dif  See  Appendix  E  for  detailed  explanation 

324,  F,l;  QENG2 

325,  Tl,  324.-QENG2  DEL 

326, +, 315, 325;  TRQENG 

327,  G, 315;  TORQ1 

328,  G, 325;  TORQ2 

329,  K;  DUMMY 


♦PAR 

PARAMETERS  : 

BIX,  PI,  P2,  P3 

14,-3000 

83,3000,4.6 

149,-3000 

302,125 

313.385.16. . 45 

314.18720. . 45877 

325.18719. . 3 
327, 3 . 18E-3 
328, 3 . 18E-3 

♦FUN 

FUNCTIONS  : 


14 


ELK  NO.  =  324 
OOOBD  PAIRS  : 


Set  Engine  2  torque  (ft-lb) 


0,18719 

1.7,18719 

1.71,0 

500,0 

BLK  NO.  = 

MODEL  COMPLETE 

•INT 

IMTB3J  PARAMS;  LOWER,  UPPER,  INTERVAL  =  0,10,0.02 

•OUT 

O/P  BLKS 

A 

O/P  PARAMS;  %  CHANGE  REQRD,  INTERVAL  *  0.0001,0.1 

•LOG3 :  ENGCTJM 

MODEL  O/P  TO  IOG3:EN3CT.MOD 

•STO 

EM3CT.MQD  NOT  CN  DSK 


OCN,  PAR,  FUN,  CR  ALL  :  A 
•LOG2 :  ENGCT^I 

MODEL  I/P  FRCM  IOG2:ENOCr.MOD 

•LOG1 ;  1992A__0 

BLOCK  O/P  TO  IOG1: 1992A.DAT 

*GOE 

I992A.DAT  NOT  ON  DSK 
•*  FUNNING  ** 

I992A.HEL  NOT  ON  DSK 

(s)  (Ft)  (Ft/fn)  (Knots) 

Time  Alt  RoC/D  Speed  Slip  %Torq 
.0  3000.  580.8  81.9  -.00  30. 

%Colctv  %Cyclic 

F-A  Lat 

-107.94  6.80  -3.06 

A.S.E.  CJiannels  ; 

PITCH  -  On 
ROLL  -  On 
YAW  -  On 
ALT  HOLD  -  Off 
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A.S.W.  Mode  : 
TRAN 


[Cyclic 

.0 

Trim  - 

3000. 

EDGE] 

580.8 

81.9 

-.00 

30 

1.0 

3010. 

580.9 

81.9 

.00 

30 

2.0 

3019. 

577.1 

81.9 

.28 

30 

3.0 

3028. 

497.0 

81.8 

2.11 

53 

4.0 

3036. 

465.6 

81.7 

-.53 

61 

5.0 

3044. 

460.3 

81.7 

-.92 

54 

6.0 

3052. 

446.9 

81.6 

-.39 

55 

7.0 

3059. 

439.0 

81.6 

-.45 

56 

8.0 

3066. 

436.1 

81.6 

-.45 

55 

9.0 

3073. 

433.5 

81.6 

-.37 

55 

10.0 

3081. 

431.9 

81.6 

-.35 

55 

RUN  CPU 

TIME  : 

22. 

11  Sec. 

*EXI 

Fortran  program  executed  STOP  statement  0 

5 .  COMPARING  FLIGHT  DATA  WITH  SEA  KING  MODEL  RESULTS 

As  well  as  being  able  to  obtain  output,  in  tabular  or  graphical  form,  of  either  the  flight 
data  or  Sea  King  model  results  separately,  the  program  TRANS  (Refs  1-3)  may  also  be 
used  to  overlay  plots  from  both  for  comparison  purposes.  An  example  is  now  given  on 
comparing  processed  flight  data  (15018.DAT)  obtained  in  Section  3  with  corresponding 
model  data  (1518A.DAT)  obtained  in  Section  4.2.1. 

Though  TRANS  allows  plot  limits  to  be  set  automatically,  this  is  inadvisable  when 
obtaining  overlay  plots,  as  the  scales  are  determined  by  the  initial  data  file  and  may  not  be 
appropriate  for  subsequent  data  files  plotted  as  overlays.  It  is  recommended  therefore  that 
the  plot  limits  be  set  using  the  SCAle  command,  either  creating  the  limits  interactively  or 
reading  them  in  from  existing  scale  files  TRANS.SCA  and  REPT.SCA.  The  file 
TRANS.SCA  is  used  by  the  first  data  file,  and  the  file  REPT.SCA,  which  should  be  a  copy 
of  TRANS.SCA,  is  used  by  subsequent  data  files  to  be  plotted  as  overlays.  For  the  case 
here,  the  file  HOVER.SCA,  which  is  on  ARL  magnetic  tape  M228,  should  be  renamed 
TRANS.SCA  and  then  copied  as  REPT.SCA. 

To  avoid  the  tedious  retyping  of  block  numbers  that  are  commonly  required  to  be  output 
in  groups,  it  is  advisable  to  use  a  TRANS. BLK  file,  which  contains  lists  of  each  group  of 
blocks  to  be  output  together  (usually  eight  variables  in  'strip'  plot  form  to  a  page).  For  the 
present  case,  the  standard  TRANS. BLK  file  (on  tape  M228),  which  was  originally 
generated  using  the  LISt  command  (a  TRANS  command  rather  than  the  Elxsi  operating 
command),  is  used. 

Given  the  files  HOVER.SCA  (renamed  to  TRANS.SCA  and  REPT.SCA)  and 
TRANS.BLK,  which  are  listed  in  Appendix  F,  TRANS  can  now  be  run  as  follows: 
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:  TRANS 

(TRANS  version  date  11 -MAR-8  6] 

I/p  FILENAME  =  15018 

15018  -  Hover  -  Pedal  Pulse  Port  -  ASE  on  -  Fit  4 


I/P  FILE  RECORDED  CN  09 -Jan- 87  AT  14:56:57 


IOTB3J  IMP  *»  .OOOOE+OO;  FUN  CPU  TIME  =  24.98  SBC. 

TIME  FRCM  .OOOOE+OO  TO  5. OOOOE+OO  IN  STEPS  CF  J..6667E-02 

‘SPA 

IS  SPACIN3  BETWEEN  PLOTS  REQRD  :  N 
*SCA 

BUK  NO.  -1  DENOTES  HBEP  VARIABLE 

ARE  PLOT  SCALE  LIMITS  TO  BE  BEAD  FBCM  DSK  :  1 
IS  TTY  LISTING  CF  LIMITS  RBQRD  :  N 
ABE  (EDIFICATIONS  RSCPD  :  N 
*PLS 

IS  GBAPHICS  OUTPUT  TO  SCREEN  REQUIRED  :  M1 
[PLS/O  Output,  for  this  run,  going  to  DSK:15018.PLT 

STRIP  PLOTS  : 

BLKS 

LI ,  L2 ,  L3 

TO  SPECIFY  NO.  CF  X  UNITS/DCH,  TYPE  0  TOR  X 

LENGTH  OF  AXES  in  INCHES;  X,  Y  =  0,1 
MIN  X,  NO.  OF  X  UNITS/INCH  =  0,1 
ARE  SYMBOLS  RBQFD  FOR  PLOTS  :  N 
LINE  KEY  <0  GIVES  DEFAULT)  -  0 
*GOE 

*»  FUNNIMS  ** 


*REP 

I/P  FILENAME  =  1518A 

ABE  SYMBOL  BEQBD  FOR  PLOTS  :  N 


1  n  should  be  noted  that  when  using  the  REPeai  command,  the  user  should  specify  that  no  graphics  output 
to  the  screen  is  required. 
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LINE  KEY  (0  GIVES  DEFAULT)  =  0 
**  PUNNING  ** 


I/P  FILENAME  = 

Fortran  program  executed  STOP  statement  1 


The  above  run  creates  the  metafiles  15018.PLT,  15018.PL0,  and  15018.PL1  (since 
three  pictures  are  requested  representing  the  groups  LI,  L2,  and  L3).  The  output  of  these 
files  obtained  on  the  Zeta8  plotter  using  the  shellfile  'Plot1,  is  shown  in  Appendix  G.  This 
shellflle,  which  is  on  ARL  magnetic  tape  M228  and  is  listed  below,  is  designed  so  that 
Zeta8  plots  of  TRANS  files  are  in  inches. 

.LIST  PLOT 

—  SHELLFILE  TO  PLOT  FILES  CREATED  BY  'TRANS '  ON  ZETAS 

—  FRAME  SIZE  ENSURES  MAXIMJM  USE  OF  PAPER  WIDTH 

—  FILENAME  AND  PICTURE  NUMBER  MAY  EITHER  BE  INCLUDED 

—  ON  OCM4AND  LINE  OR  WILL  BE  PROMPTED  FOR 
parm  file 

pam  no 
label  get_file 

—  If  filename  has  been  specified  go  to  label  getjpicture_no 

—  otherwise  prompt  for  filename 
if  [file]  then 

else 

declare  file 

echo  'Filename  :  '  -n 

—  n*  stops  cursor  on  same  line  as  echo 

—  search  line  for  1  string  and  place  it  in  file 
set  file  [find  ?*  occurrenoe=l] 

end  if 

—  if  file  exists  go  to  label  go t_pic t ure_no ,  otherwise  prompt 

—  for  another  filename 

if  [FileExists  [file]]  then 
else 

echo  'Filename  '[file]'  does  not  exist' 
forget  file 
goto  get_file 
end  if 

label  getj?icture__no 

—  if  picture  ranter  is  specified,  plot  file,  otherwise  prompt 

—  for  picture  ranter 
if  [no]  then 

else 

declare  no 

echo  'Picture  ranter  :  ’  -n 
set  no  [find  ?*  occurrence^!] 
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end  if 

—  if  no  is  not  a  nunfoer,  then  ptcnpt  for  picture  nuiter 

if  (echo  [no]  I  find  %0<9*$]  then 

else 

echo  ' Invalid  picture  nunber  1 [no] 
forget  no 

goto  get_picture_no 
end  if 

plot.zt8  [file]  picno=[no]  frame=534, 267 
printers 

Plot  operates  as  shown  below.  Note  that  Plot  will  prompt  for  the  filename  and/or 
picture  number,  if  they  are  not  included  in  the  command  line. 

:Plot  15018. PLT  1 


jobID 

state 

user 

title 

size 

pri 

time 

destination 

89 

active 

ae  ,arrtey 

15018  - PLT 

436 

norm 

15:22 

cc, zeta8 

6.  CONCLUDING  REMARKS 

Perhaps  the  most  significant  achievement  in  developing  the  Sea  King  model  has  been 
the  extent  to  which  it  has  been  validated  against  flight  data.  The  capabilities  of  the  model, 
as  well  as  limitations,  are  consequently  reasonably  well  understood.  However,  this  and 
any  further  validation  requires  a  knowledge  of  how  to  process  the  ARL  flight  data,  as  well 
as  how  to  run  the  model  using  these  data  as  inputs  in  order  to  ensure  compatibility  between 
theoretical  and  experimental  results.  To  keep  alive  this  capability,  especially  following 
transfer  of  the  model  and  flight  data  from  the  obsolete  DEC  System  10  to  the  Elxsi  6400,  it 
has  been  necessary  therefore  to  fully  document  all  aspects  of  the  model  in  the  form  of  these 
user  manuals.  Inclusion  here  of  a  comprehensive  catalogue  of  the  computer  records  of  the 
flight  data  has  been  considered  necessary  in  view  of  current  expectations  on  the  value  of  the 
flight  data  in  future  studies,  both  in  the  application  of  parameter  identification  techniques  to 
helicopter  modelling  and  in  the  validation  of  a  Sea  King  analysis,  rather  than  simulation, 
model  being  developed  using  the  CAMRAD  (Comprehensive  Analytical  Model  of 
Rotorcraft  Aerodynamics  and  Dynamics)  code. 
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APPENDIX  E 

Model  Modifications  for  Single  Engine  Cut 

A  description  of  the  standard  engine  model  as  simulated  in  CSMP-10(ARL)  is  given  in 
Fig.  7  of  Ref.  6.  The  changes  to  this  model  (via  CONfiguration,  PARameter,  and 
FUNction  commands)  are  shown  in  Figure  E.  1  of  this  appendix. 

In  order  to  model  the  response  of  a  single  engine,  the  model  is  modified  so  that  the 
response  of  the  cut  engine  (Engine  2)  is  fed  in  as  a  FUNction  block,  with  a  first  order  lag 
simulating  the  torque  decay. 

Referring  to  Figure  E.  1 ,  the  detailed  changes  are  as  follows: 


Blk 

302  Gain  of  250  is  halved  to  125 

313  Initial  condition  halved;  time  constant  unchanged  from  standard  value 

314  Initial  condition  halved;  time  constant  unchanged  from  standard  value 

325  Initial  condition  is  half  of  total  engine  torque  (lb  ft);  time  constant  =  0.3  s 
(found  from  inspection  of  flight  data) 

327  Gain  of  3.18  x  10  3  to  convert  torque  (lb  ft)  to  %*  (see  below) 

328  Gain  of  3.18  x  10‘3  to  convert  torque  (lb  ft)  to  %  (see  below) 


Calculation  of  Gains  for  Blks  327,  328: 

From  Ref.  20, 

111%  torque  =  2778  SHP  at  103%  Nr 
and  100%  Nr  =  203  rpm 

103%  Nr  a  1.03  x  203  =  209  rpm  =  21.89  rad/s 

=>  111%  torque  =  (  27'Ji  g9  ~  )  *  698 10  lb  ft  for  two  engines 
Hence,  for  one  engine, 

Gain  =  (  0.5  x  69810  )  =  318  x  10-3 


1  %  is  unit  displayed  on  cockpit  gauge  and  is  based  on  a  nominal  value  (used  for  torque  and  rotor  revs. 
Nr). 
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APPENDIX  F 

Files  Required  to  Run  TRANS 


a)  Scale  File  HOVER.SCA 


Blk  No. 

Lower  Limit 

tapper  Limit 

-1 

O.OOOOE+OO 

5.0000E+00 

276 

-1.0000E+01 

1.  OOOOE+Ol 

277 

-1 . OOOOE+Ol 

1. OOOOE+Ol 

72 

-1.0000E-01 

1.0000E-01 

278 

O.OOOOE+OO 

2. OOOOE+Ol 

292 

-1. OOOOE+Ol 

1. OOOOE+Ol 

293 

-1. 0000E+01 

l.OOOOE+Ol 

295 

O.OOOOE+OO 

2.  OOOOE+Ol 

296 

O.OOOOE+OO 

l.OOOOE+Ol 

282 

-2.  OOOOE+Ol 

2. OOOOE+Ol 

283 

-2. OOOOE+Ol 

2. OOOOE+Ol 

284 

-2 . OOOOE+Ol 

2. OOOOE+Ol 

279 

-1. 0000E+01 

l.OOOOE+Ol 

280 

-1 .  OOOOE+Ol 

l.OOOOE+Ol 

281 

-1.0000E+01 

l.OOOOE+Ol 

287 

-2.  OOOOE+Ol 

2. OOOOE+Ol 

288 

-2.  OOOOE+Ol 

2. OOOOE+Ol 

16 

-4. 0000E+01 

-2.  OOOOE+Ol 

17 

-l.OOOOE+Ol 

l.OOOOE+Ol 

18 

-1.0000E+01 

l.OOOOE+Ol 

212 

6. OOOOE+Ol 

1.0000E+02 

290 

1.0200E+02 

1 . 0600E+02 

Note:  Headings  shown  are  not  stored  in  file. 
Blk  no.  -1  represents  time. 

b)  Block  List  File  TRANS. BLK 


27  6 

L2 

277 

72 

278 

2  92 

293 

295 

296 

282 

L3 

283 

284 

279 

280 

281 

287 

288 

16 

17 

18 

212 

29 
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